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> ; ABSTRACT 

a^ . 

<^ . We report on long-term monitoring of anomalous X-ray pulsars (AXPs) using the Rossi X-ray Timing 
I Explorer (RXTE). Using phase-coherent timing, we find a wide variety of behaviors among the sources, 
' ranging from high stability (in IE 2259.1+586 in quiescence and 4U 0142+61), to instabilities so severe that 
phase-coherent timing is not possible (in IE 1048.1-5937). We note a correlation in which timing stability 
in AXPs decreases with increasing z>. The timing stability of soft gamma repeaters (SGRs) in quiescence 
Q-i is consistent with this trend, which is similar to one seen in radio pulsars. We find no significant pulse 
Q ' morphology variations in any AXP in quiescence. We considered high signal-to-noise average pulse profiles 
for each AXP as a function of energy. We show that, as in the timing properties, there is a variety of different 
behaviors for the energy dependence. We also used the monitoring and archival data to obtain pulsed flux 
time series for each source. We have found no large changes in pulsed flux for any source in quiescence, 
. and have set la upper limits on variations ~20-30% depending on the source. We have recently discovered 
K*N ' bursts from the direction of two AXPs: IE 1048.1-5937 the most SGR-like AXP, and IE 2259.1+586 the 
d most rotationally stable AXP. We compare the temporal, spectral and flux properties of these events to those 
of SGR bursts, and show that the two phenomena are very similar. These results imply a close relationship 
between AXPs and SGRs, with both being magnetars. 
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INTRODUCTION 

The nature of anomalous X-ray p ulsars (AXPs) was a mystery for the past 20 years, since the discovery 
of the first example (IE 2259.1+586; iFahlman and Gregorv. 198ll ). There are currently only five confirmed 
AXPs, all of which are located in the Galactic plane, with two at the geometric centers of supernova remnants 
dPahlman and Gregorv. 198llVasisht and Gotthelf. 1997 ) . AXP characteristics can be summarized as follows 
fsee llsrael et al. 20021 . for a review): they exhibit X-ray pulsations in the range ~6~12 s; they have pulsed 
X-ray luminosities in the range ~ 10^^ — 10^^ erg s~^; they spin down regularly; they have spectra that 
are characterized by therma l emission of kT ^0.4 keV plus a hard power-law tail. They have faint or 
no optical/IR counterparts ( Hulleman et al.. 2000slHulleman et al.. 2000lJlsrael et al.. 2002h . The optical 
counterpart of one AXP, 4U 0142+61 is reportedly pulsed (|Kern and Martin. 20021 ). 

AXPs are called "anomalous" because it has been unclear what powers their radiation. They are not 
rotation-powered as their observed X-ray luminosities are much greater than the rate of loss of rotational 
kinetic ener gy inferred from their spin-do wn. They were long thought to be accreting from a low mass com- 
panion (e.g. Mereghetti and Stella. 19951 ). However this model is difficult to reconcile with observations: the 
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absence of Doppler shifts even on short time scales (e.g. iMereghetti et al., 19981) , the absence of a detec table 
optical/IR companion or bright accretion disk fsee iMereehetti and Stella. 1995Hulleman et al.. 2000al ). the 
apparent associations with supernova remnants, that AXP spectra are very different from those of known ac- 
creting s ources, and that is gen erally s maller than in known accreting sources, all are inconsistent with this 



scenario. 



Chatteriee et a,]. (200(1 ) (see a.lso lva.n Paradiis et al.. 199,iCorbet et al.. 199,'ih proposed that AXPs 



are accreting from a fossil disk made of fall-bac k material from the supernova. However, this model overpre- 

di ets the observed optical/IR flux from the disk ( Perna et al.. 200(iHulleman et al.. 2000slHulleman et al.. 2000lJKaplan e 



Thompson and Duncan (199fil ) suggested that the main source of free energy for AXP emission is from 



the magnetic field itself. In this model, AXPs are young, isolated, highly magnetized neutron stars or 
"magnetars". Prior to our RXTE monitoring campaign, the primary evidence in favor of this model for 
AXPs was the inferred strength of the surface dipolar magnetic field required to slow the pulsar down in 
vacuo, assuming magnetic dipole braking as in radio pulsars. These fields are in the range ~ 10^^ — 10^^ G. 
The spin-down ages in this model, inferred assuming a small birth spin period, are in the range ~8-200 kyr. 
This suggested youth is further evidenced by the associations with supernova remnants, and from the location 
of AXPs in the Galactic plane. 

The identification of AXPs with magnetars was further strongly motivated by the similarity of AXP 
emission to that of the soft gam ma-rav repeaters (SGRs) in quiescence. Specifical ly, the latter have similar 
pulse periods, are spinning down ( Kouveliotou et al.. 199dKouveliotou et al.. 19991 ). and have X-ray spect 
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that are comparable to, though somewhat h arder than, those of the AXPs, at least when not in outburst 
( Mereghetti et al.. 2nnfiKulkarni et al.. 200'^ . independent evidence for the ultra-high magnetic fields exists 



in SGRs; for exar nple, a ^ 10^^ G magnetic fiel d is required to confine the radiation that is seen following 
major outbursts ( Thompson and Duncan. 1995h . Indeed the only strong distinction between the emission 



from SGRs and AXPs appeared to be the fact that SGRs exhibit occasional and repeating short bursts, 
while AXPs did not. 

RESULTS 

The results presented here were obtained using the Proportional Counter Array fPCA: ljahoda et al.. 199(} ) 
on board NASA's Rossi X-ray Timing Explorer (RXTE). The PCA consists of an array of five collimated 
xenon/methane multi-anode proportional counter units operating in the 2-60 keV range, with a total 
effective area of approximately 6500 cm^ and a field of view of ~ 1° FWHM. A program to monitor AXPs 
regularly using RXTE was begun in 1996 in order to determine their long-term timing, pulsed fiux , and 
pulse profile stabilities ( Kaspi et al.. 199tiKaspi et al.. 200ClKaspi et al.. 200]lGavriil and Kaspi. 2002h . As 



part of this prograrn , motivated by the existence of SGR bursts, we also searched the AXP data for bursts 
( Gavriil et al.. 20021 ). Our observations consist primarily of short snapshots taken on a monthly basis. In 



addition, we used a handful of archival observations; the exposures in these observations vary. We used 
the GoodXenonwithPropane data mode, which records photon arrival times with l-/xs resolution and bins 
ener gies into one of 256 channels. Below, w e sum marize the res ults of this monitoring campaign thus far 
(see Kaspi et al.. 200]lGavrnl and Kaspi. 2002.Gavriil et al.. 20021 . for details). 



Phase- Coherent Timing 

In the timing analysis, each binned time series was epoch-folded using the best estimate frequency de- 
termined initially from either a periodogram or Fourier transform (though later folding was done using the 
timing ephemeris determined by maintaining phase coherence; see below). Resulting pulse profiles were 
cross-correlated in the Fourier domain with a high signal-to-noise template created by adding phase-aligned 
profiles from previous observations. The cross-correlation returns an average pulse time-of-arrival (TOA) 
for each observation corresponding to a fixed pulse phase. The pulse phase (p at any time t can be expressed 
as a Taylor expansion, 

(j){t) = ^{to) + uoit - to) + - to)^ + - to)^ + ■■■, (1) 

2 D 

where = 1/P is the pulse frequency, z> = dv/dt, etc. and subscript '0' denotes a parameter evaluated at the 
reference epoch t = to- The TOAs were fit to the above polynomial using the pulsar timing software package 
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Fig. 1. Arrival time residuals for 4U 0142+61 with u 
and z> subtracted. The RMS residual here is < 0.3% of 
the period, indicating great rotational stability. How- 
ever, we realize that the reduced is much greater 
than unity, this is under investigation. 
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Fig. 2. Spin history for IE 1048.1-5937. The points 
represent past measurements of the frequency of the 
pulsar (see Oosterbroek et al. 1998 and references 
therein, Paul etal. 2000 and Baykaletal. 2000). Solid 
lines represent our RXTE observations (Kaspi et al. 
2001). Panel A shows the observed frequencies over 
time. The dotted line is the extrapolation of a and 
j> measured from a phase-connected interval in 1999. 
Panel B shows the difference between the ephemeris 
indicated by the dotted line and the data points. 



TEMPO^. Unambiguous pulse numbering is made possible by obtaining monitoring observations spaced so that 
the best-fit model parameters have a small enough uncertainty to allow prediction of the phase of the next 
observation to within ~ 0.2. Typically this requires two closely spaced observations (within a few of hours 
of each other) followed by one spaced a few days later, and regular monitoring thereafter, as long as phase 
coherence can be maintained. 

Using phase coherent timing we have shown that AXPs can be phase connected over long intervals, 
implying they can spin down with impress ive stability, which is at odds with what is seen in most X -ray 
pulsars that definitely have accretion disks ( Kaspi et al.. IDQ^Kaspi et al.. 200ClGavriil and Kaspi. 20021 . see 



Figure However t his is not true of al l our sources. IE 1048.1-5937 shows large deviations from a simple 



spin down (Figure 121 iKaspi et al.. 2001 ). However, the very noisy timing behavior seen in IE 1048.1-5937 



shows no evidence for correlated flux variations, as are expected in fall-back disk models. 
Pulse Morphology 

Motivated by a p revious report of a signiflcant pulse morphology change in an AXP (IE 2259.1-1-586; 
Iwasawa et al.. 19921 ) we used our RXTE data to search for more such events. Morphology changes are 



expected in the magnetar model during outbursts or torque variations. 

To search for pulse profile changes, the profiles were first phase aligned using the templates and the same 
cross-correlation procedure used for timing. Each data profile was fit to a high signal-to-noise template by 
adjusting the amplitude and the offset in order to minimize a statistic. The resulting data profile was 
subtracted from the template to yield "profile residuals" for each observation. The procedure was repeated 
for each pulsar with different binning in order to have sensitivity to a variety of types of pulse profile changes. 



^http://pulsar. Princeton. edu/tempo 
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We have not detected any large pulse profile variations, except for IE 2259.1+586 during a major outburst 
(see below). This justifies our other analysis procedures which assume a fixed profile. However, in a handful 
of observations, we have found low-level pulse profile changes, at the ~ 3a level for all sources. 

It is difficult to set quantitative upper limits on the amplitude of pulse profile changes to which we were 
sensitive, as these depend on the shape of the change, and vary depending on the length of the observation. 
Typically, RMS profile residuals are ~20% of the pulse peak, although this varied from 4-40%. 

Pulsed Flux and Spectra 

We have also used our RXTE observations to monitor the pulsed fiux of each AXP. By considering the 
stability of pulsed fiuxes, we test the two competing AXP models. In accreting systems, torque fluctuations 
are generally accompanied by variation in the mass accretion rate and hence luminosity, and vice versa. In 
the magnetar model, abrupt flux changes could occur in analogy with SGR-like outbursts. Given the large 
field-of-view of the PCA, the low count rates for the sources relative to the background, and the fact that, 
for example, IE 2259.1+586 is in a supernova remnant, total fiux measurements are difficult with our RXTE 
data. Instead, we have determined the pulsed component of the fiux, by using the off-pulse emission as a 
background estimator. 

Data from each observing epoch were folded at the expected pulse period as was done for the timing 
analysis. However, for the flux analysis, 16 phase bins were used across the pulse. For each phase bin, 
we maintained a spectral resolution of 128 bins over the PCA range. Given the broad morphologies of 
the average pulse profile, only one phase bin could be used as a background estimator. The pulse profiles 
were phase aligned, so that the same off-pulse bin was used for background in every case. The remaining 
phase bins were summed, and their spectral bins regrouped using the FTODL grppha. The regrouped, phase- 
summed data sets, along with the background measurement, were used as input to the X-ray spectral 
fitting software package XSPEC^. Response matrices were created using the FTOOLs xtefilt and pcarsp. 
Because of the limited statistics, fitting a two-component model was not practical, so we used a simple 
photoelectrically absorbed power law, at first holding only Nh fixed. For all sources we found that the 
photon index T was constant within the uncertainties; we therefore held it fixed at its mean value. To 
extract a pulsed fiux at each observing epoch we refit each spectrum by varying only the normalization. 
Other than after a major outburst from IE 2259.1+586 (see below), we found no evidence for any large 
variability in the pulsed flux within 20%-30% consistent with the systematic errors. AXP spectra are 
generally best flt by a two-component model consisting of a photoelectrically absorbed blackbody with 
a hard power-l aw tail (Isra el et al . 1999). Whether these two components are physically distinct is an 
open question ( Ozel et al.. 20fll,Thompson et al.. 20021 ). To investigate this, we compared the pulse proflle 



morphology of the AXPs in two energy bands. Thus far, we have found that pulse proflle energy-dependence 
is generally subtle, but varies from so urce to source. We found no clear patterns when considering average 
pulse profiles as a function of energy dGavriil and Kasui. 20oi) . 



SGR-like X-ray Bursts from IE 1048.1-5937 

In the magnetar model for SGRs ( Thompson and Duncan. 19961 ). bursts are a result of sudden crustal 



yields due to stress from the outward diffusion of the huge inter nal magnetic field. Such yields cause crust 
shears which twist the external magnetic field, releasing energy. Thompson and Duncan (19961 ) who, upon 



suggesting that AXPs are also magnetars, predicted X-ray bursts should eventually be seen from them. In 
order to confirm this pr ediction we searc hed our RXTE data for bursts. Our burst searching algorithm is 



summarized below (see iGavriil et al.. 2002 . for details). 



Time series were created with 31.25-ms resolution from photons having energies in the range 2-20 keV 
for each PCA Proportional Counter Unit (PCU) separately, using all xenon layers. Photon arrival times 
at each epoch were adjusted to the solar system barycentre. The resulting time series were searched for 
significant excursions from the mean count rate by comparing each time bin value with a windowed 7-s 
running mean. Bursts were identified assuming Poissonian statistics, and by combining probabilities from 
the separate PCUs. 

We discovered two short, highly significant X-ray bursts, separated by 16 days, from the direction of 



http://xspec.gsfc.nasa.gov 
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Fig. 3. Left Panel: Background subtracted 2-20 keV light curves for the first burst displayed with 31.25-ms time 
resolution. Right Panel: Same but for second burst. 



the AXP IE 1048.1-5937 in data from Fall 2001 (|Gavrii1 et a,1.. 2002l V The first occurred during a 3-ks 
PCA observation obtained on 2001 October 29, and the second during a 3-ks observation obtained on 2001 
November 14. We searched a total of 442 ks of PCA observing time from 1996-2002, and no other significant 
bursts were found toward IE 1048—5937. The burst profiles are shown in Figure |31 

We fit the spectra from the first 1 s of each burst to two one-component models, a power law and a black 
body. Continuum models provided an adequate characterization of the Burst 2 spectrum but not of the 
Burst 1 spectrum. As seen in Figure|lJ the spectrum for the 1 s after the Burst 1 onset exhibits a feature near 
14 keV. This feature is clear in all binning schemes and is prominent throughout the first ~1 s of the burst. 
No known PCA instrumental effect produces a feature at this energy (K. Jahoda, personal communication). 
Intriguingly, in some spectral binning schemes, we saw hints of features in the spectrum consistent with 
additional multiples of 7 keV, suggesting cyclotron lines. However these were not statistically significant. 

Due to the wide (~1°) field-of-view (FOV) and lack of imaging capabilities of the PCA, we cannot verify 
that the bursts originated from the location of the AXP. The low peak X-ray fluxes of the events (see Table 
1) preclude determining the source's location using data from other, better imaging instruments that were 
contemporaneously observing the X-ray sky, such as the RXTE All Sky Monitor, or the Wide Field Camera 
aboard BeppoSAX. We must therefore consider other possible origins from the bursts before concluding they 
were from the AXP. 

The bursts' short rise times (Table 1) require emission regions of less than a few thousand km, implying 
a compact object origin. Type I bursts from an LMXB in the same FOV as IE 1048.1-5937 are unlikely to 
explain our observed bursts because (i) the burst rise times are much shorter than those of Type I bursts; (ii) 
the burst spectra are much harder than those of Type I bursts; (iii) Burst 2 shows no evidence for spectral 
softening with time and no Type I burst has ever exhibited a spectral feature like the one detected in Burst 
1; (iv) the bursts are extremely faint, implying a source locati on well outside the Milky Way for Type I 
burst luminosities (v) there are no known LMXBs in the FOV (|Lhi et al . 200ll 'l. Furthermore, the bursts 
we have observed are unlikely to be Type II bursts from an unknown X-ray binary in the PCA FOV because 
(i) of the rarity of such events (observed thus far only in two sources, both accreting binaries); (ii) Type II 
bursts have longer rise times than do our bursts; (iii) no Type II burst has exhibited a spectral feature like 
that seen in Burst 1. Classical gamma-ray bursts (GRBs) sometimes exhibi t prompt X-rav em ission that 
can have temporal and spectral signatures similar to those we have observed (|Heise et al.. 20011 ). However, 
the likelihood of two GRBs occurring within 1° of each other is small, and GRBs are not known to repeat . 

The observed burst properties are in many ways similar to those seen from SGRs (|Gogiis et al.. 200lh . 
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The fast rise and slow decay profiles are consistent with SGR time histories, as are the burst durations. 
Both AXP and SGR bursts are spectrally much harder than is their quiescent pulsed emission. The burst 
peak fluxes and fluences fall within the range seen for SGRs, and the spectrum of Burst 2 is consistent 
with SGR burst spectra of comparable fluence. Burst 1 has characteristics unlike nearly all SGR bursts, 
specifically i ts long tail fTable 1) and spectral feature. However, we note that one event from SGR 1900+14 
was shown ( Ibrahim et al.. 200]IStrohmaver and Ibrahim. 200nl ) to possess each of these properties. The 



marginal increase in the pulsed flux that we observed at the bu rst epochs is consistent with SGR pulsed 
flux increases seen during bursting episodes ( Woods et al.. 200lh . Finally, the fact that in spite of several 



years of monitoring, the only two bursts detected occurred within two weeks of each other suggests episodic 
bursting activity, the hallmark of SGRs. Thus, the characteristics of these events match the burst properties 
of SGRs far better than any other known burst phenomenon. 

The large 14-keV spectral line in the first burst in intriguing. An electron cyclotron feature at this energy 
E implies a magnetic field of B = 2TrmcE/he ~ 1.2 x 10^^ G (where m is the electron mass, c is the 
speed of light, h is Planck's constant, and e is the electron charge), while a proton cyclotron feature implies 
B ~ 2.4 X 10^^ G. The former is significantly lower than is implied from the source's spin-down and is typical 
of conventional young neutron stars, rather than magnetars. The latter is higher than is implied by the 
spin-down yet reasonable for the magnetar model as the spin-down torque is sensitive only to the dipolar 
component of the magnetic field. 

A Major Outburst from IE 2259.1-|-586 

On 2002 June 18, fortuitously during one of our RX TE monitorm^ observations of IE 2259.1+586, over 
80 bursts were detected in the span of 14 ks ( Kaspi and Gavriil. 200jKasni et aJ.. 2002l V Figure 13 shows the 



0.125 s light curve. IE 2259.1+586, in contrast to IE 1048.1-5937, had been an extremely stable rotator 
and had shown no signs of "volatility." The association of the bursts toward IE 1048.1-5937 with the AXP 
could not be unambiguously verified, however, the bursts seen here (Figure [51 top panel) clearly originated 
from IE 2259.1+586. During and following the outburst, the pulsed flux increased (see Figure ISJ and there 
was a significant pulsed morphology change, such that the relative amplitude of the leading and trailing 
pulse were swapped (see Figure IH)). 

Target-of-opportunity observations were made by many major observatories in addition to RXTE at the 
time of this outburst. No radio emission was detected in spite of sensitive VLA observations, nor were 
any bursts seen in TOO XMM-Newton data obtained roughly a week afterward. Near- infrared observations 
obtained with the Gemini North 8-m reflector on June 21 showed that the proposed IR counterpart had 
brightened by a factor of 3.4 (Kaspi et al. 2002). A few days later, it faded back to its pre-burst value. This 
represents the firs t detection of transient optical/IR emission from an isolated neutron star. Subsequently, 
Israel et al. (200^ reported possible IR variability from the proposed counterpart to IE 1048.1—5937, which 



they suggested might be related to its bursting behaviour. The origin of IR emission is currently not 
addressed in the magnetar model. 

Pulse morphology changes can occur, according to the magnetar model, during outbursts as in the SGRs, 
as this is when magnetic reconfigurations occur. Indeed during the outburst of IE 2259.1+586, significant 
pulse morphology chan ges were seen, such t hat the relative amplitude of the double peaked profile was 
swapped (see Figure 1^ . Ilwasawa et al. ('1992h also reported a significant change in the pulse morphology of 
IE 2259.1+586 in GINGA observations obtained in 1990, such that the leading pulse had amplitude roughly 
half of the trailing pulse. In both instances the pulse profile returned back to its long-term state. 

Pulsed flux enhancements in both the pulsed and persistent fiux are observed in the SGRs during bursting 
episodes and are thought to be a result of back heating of the surface and outer crust by the burst emission in 
the magnetosphere. The pulsed fiux of IE 2259.1+586 increased by a factor of ~ 4 during its bursting episode 
(see Figureini). We noted a marginal (^ 3cr) increase in the pulsed fiux from IE 1 048.1—5937 that com menced 
with the observation in which Burst 1 was detected, and which lasted ~4 weeks ( Gavriil et al.. 20021 ). There 



have been past clair ns of X-rav fiux c hanges, bv up to a factor of ~10 in two AXPs (IE 2259.1+586 



and IE 1048.1-5937; iBavkal and Swank . 199(: Oosterbroek et al.. 19 98). although the reported fiuxes were 



measured with different instruments having very different spectral responses and angular resolutions. If 
correct, these past claims of increased fiux might imply previous bursting activity. 
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Parameter 


Value 




Burst 1 


Burst 2 


Burst day, (MJD) 


52211 


52227 


Burst start time, 


0.2301949(24) 


0.836323379(68) 


(fraction of day, UT) 






Burst rise time, tr (ms) 


2lt^ 


5.912;° 


Burst duration, Tqq (s) 


t-1 +28 
5III9 


n n+4 9 


igo iluence (counts j 


/I Q _L 1 1 Q 

4oo ± iio 


lUi ± io 


1 90 nuence i^xiu eig cm j 


on 01 AO 

20.3 ± 4.8 


5.3 ± 1.2 


Peak flux for 64 ms 


01+45 
31-10 


26lf 


f in — 10 —1 

(x 10 erg s cm j 




Peak nux for 




114I23 


(xlO ^ erg s cm ) 






power law index 


o.89ii:?i 


1 Qo+u.y5 

-'-•"J°-0.62 


power law flux 


^•'-'-1.8 


^■U-0.8 


(xlO"^" erg s~^ cm~^) 




line energy (keV) 


13.9 ±0.9 




line width, a (keV) 


9 9+1.3 
^•^-1.0 




line flux 


q+2.2 




(xlO~^° erg s~^ cm~^) 




reduced (i^) 


1.24 (15) 


0.77 (5) 



Table 1. AXP IE 1048.1-5937 burst properties. The burst rise tr times 
were determined by a maximum likelihood fit to the unbinned data using 
a piecewise function having a linear rise and exponential decay. The burst 
duration, T90, is the interval between when 5% and 95% of the total 
2-20 keV burst fluence was received. All fluences and fluxes are in the 
2-20 keV range. T90 fluences in cgs units were calculated assuming a 
power-law spectral model. For all spectral fits, nu was held fixed at 
1.2 X 10^2 cm~2, th e value determined from recent XMM observations 
( TiengQ et al.. 20o3 ). Spectral modeling was done using photons in the 
2-40 keV range. Uncertainties are 68% confidence intervals, except for 
those reported for the cgs-unit fluences and fluxes, as well as the spectral 
model parameters, for which we report 90%) confidence intervals. See 
Gavriil et a for details. 




10 20 
channel energy (keV) 



Fig. 4. X-ray spectrum in the 
2-40 keV range for the 1 s after 
the onset of Burst 1. The best fit 
power-law plus Gaussian line model 
is shown as a solid line. The F-test 
shows that the addition of a line of 
arbitrary energy, width and normal- 
ization to a simple power law model 
improves the fit significantly, with a 
chance probability of this occuring of 
0.0032. Monte Carlo simulations in 
XSPEC were done to verify this con- 
clusion. 



DISCUSSION 

We have presented a variety of observational results for all AXPs using regular monitoring observations and 
archival observations from the PCA aboard RXTE. There is now strong evidence that AXPs are magnetars. 
The great rotationally stability in some AXPs is at odds with an accretion scenario, as is the absence of 
pulsed flux or pulse morphology changes correlated with torque. Extended periods of spin-up would only, 
though not necessarily, be expected in an accretion scenario, however all AXPs have exhibited relatively 
steady spin-down. Furthermore, and perhaps most compellingly, accretion models cannot account for the 
bursting episodes seen in two AXPs. Interestingly, there is a correlation between timing stability and v in 
our data. The sources with by far the smallest i>'s, are the most stable in quiescence (at least during ou r 
observations), while those with large z>'s are less so. The SGRs, being even less stable ( Woods et al.. 20o3 ). 
have even larger z>'s, in agreement with this trend. If correct, the trend suggests a continuum of timing 
properties between the AXP and SGR populations, lending additional support to the connection between 
them. 

The identification with magnetars is further strongly motivated by the similarity o f the AXP emission to 

that of the SGRs. Specifically, the latter have similar pulse periods, are spinning down ( Kouveliotou et al.. 199feouveliot" 
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Fig. 5. The time evolution of IE 2259+586 during 
its burst active episode on 2002 June 18 as observed 
with the RXTE PCA. Top: lightcurve in 2-20 keV X- 
rays with 0.125 s time resolution. Bottom: The pulsed 
2-10 keV flux. The crosses represent the pulsed flux 
in the 2-20 keV band and the diamonds represent the 
unabsorbed persistent flux. 
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Fig. 6. Average pulse profiles of IE 2259.1-1-586 in the 
2.5-9.0 keV band. Two cycles are plotted for clarity. 
(A) Average profile before the outburst (total exposure 
time: 764 ks). (B) Average profile during the outburst, 
with bursts omitted (total exposure time 11 ks). Dur- 
ing the outburst, the amplitudes of the peaks relative 
to the pre-outburst profile are clearly reversed. The 
relative phase displayed is that successfully used in our 
timing analysis. 



and have X-ray spectra that are comparable to, though somewhat harder than, those of the AXPs, at least 
when not in outburst. With the detection of X-ray bursts from two AXPs the two classes are now simi- 
lar in that respect as well. These very rec ent results appear to have seal ed the connection between SGRs 
and AXPs, as was originally predicted by Thompson and Duncan (199()l ). However, the identification of 



the two classes objects as magnetars has left some open questions. For instance, w hat if anvthing distin - 
guishes AXPs and SGRs? Are AXPs the progenitors of SGRs or vice versa? Recentlv lKulkarni et al. 



have shown that the X-ray spectrum of the quiescent counterpart to SGR 0526-66 has a power-law index 
—3, similar to those seen in AXPs, but softer than those of SGRs. As SGR 0526-66 has been inactive 
for almost two decades, this suggests SGRs and AXPs are one and the same object, and may have even 
cycle back and forth between each flavor. Another noteworthy question is what is the magnetar/radio 
pulsar conne ction? Recent obse rvations of radio pulsars having apparent surface dipolar magnetic fields 
comparable ( Camilo et al.. 2000I ). and in one very new case (McLaughlin et al., private communication). 



surpassing those of some AXPs, leaves open the questi on of what phv sicallv distinguishes a magnetar from a 
rotation-powered pulsar of high magnetic field? Pivovaroff et al. (2000ii) suggested that the sources classified 
as magnetars might have much higher fields, in the form of multipole moments. Perhaps comparing the 
thermal X-ray emission from radio pulsars having a variety of magnetic fields may illuminate the situation. 
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